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Abstract: The semi-phenomeno logical improved holographic model for QCD is 
confronted with data of the pure glue, large-iV c gauge theory. After fitting two 
phenomenological parameters in the potential, the model can reproduce in detail 
all thermodynamic functions at finite temperature. It also reproduces in detail all 
known spin-0 and spin-2 glueball observables at zero temperature and predicts the 
rest of the ++ and 2 ++ towers. A similar two parameter fit in the CP-odd sector 
postdicts the correct second + ~ glueball mass, and predicts the rest of the + ~ 
tower. 
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1. Introduction 

The experimental efforts at RHIC, have provided a novel window in the physics of 
the strong interactions. The consensus on the existing data is that shortly after the 
collision, a ball of quark-gluon plasma (QGP) forms that is at thermal equilibrium, 
and subsequently expands until its temperature falls below the QCD transition (or 
crossover) where it finally hadronizes. Relativistic hydrodynamics describes very well 
the QGP 0, with a shear- viscosity to entropy density ratio close to that of M = 4 
SYM, 0. The QGP is at strong coupling, and it necessitates a treatment beyond 
perturbative QCD approaches, 0. Moreover, although the shear viscosity from M = 
4 seems to be close to that "measured" by experiment, lattice data indicate that in 
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the relevant RHIC range 1< £ < 3 the QGP seems not to be a fully conformal fluid. 

J- c 

Therefore the bulk viscosity may play a role near the phase transition || |J. The 
lattice techniques have been successfully used to study the thermal behavior of QCD, 
however they are not easily extended to the computation of hydrodynamic quantities. 
They can be used however, together with parametrizations of correlators in order to 
pin down parameters |J. On the other hand, approaches based on holography have 
the potential to address directly the real-time strong coupling physics relevant for 
experiment. 

In the bottom-up holographic model of AdS/QCD [J/J, the bulk viscosity is zero 
as conformal invariance is essentially not broken (the stress tensor is traceless). In the 
soft-wall model ||, no reliable calculation can be done for glue correlators and there- 
fore transport coefficients are ill-defined. Similar remarks hold for other phenomeno- 
logically interesting observables as the drag force and the jet quenching parameter 
[@, 0,0,0. 

Top-down holographic models of QCD displaying all relevant features of the the- 
ory have been difficult to obtain. Bottom-up models based on AdS slices [|IB[] have 
given some insights mostly in the meson sector, but necessarily lack many impor- 
tant holographic features of QCD. A hybrid approach has been advocated |14, 15 



combining features of bottom-up and top-down models. Such an approach is essen- 
tially a five-dimensional dilaton-gravity system with a non-trivial dilaton potential. 
Flavor can be eventually added in the form of Nf space-time filling DA — D4 brane 



pairs, supporting U(Nj) L x U(Nf) R gauge fields and a bi-fundamental scalar [|R 
The UV asymptotics of the potential are fixed by QCD perturbation theory, while 
the IR asymptotics of the potential can be fixed by confinement and linear glueball 



asymptotics. An analysis of the finite temperature behavior [[17], [18|] has shown that 
the phase structure is exactly what one would expect from YM. A potential with a 
single free parameter tuned to match the zero temperature glueball spectrum was 
able to agree with the thermodynamic behavior of glue to a good degree, |17]]. Similar 
results, but with somewhat different potentials were also obtained in |19|, |20| 

In |17], [18|] it was shown that Einstein-dilaton gravity with a strictly monotonic 
dilaton potential that grows sufficiently fast, generically shares the same phase struc- 
ture and thermodynamics of finite-temperature pure Yang-Mills theory at large N c . 
There is a deconfinement phase transition (dual to a Hawking-Page phase transi- 
tion between a black- hole and thermal gas background on the gravity side), which is 
generically first order. The latent heat scales as N%. In the deconfined gluon-plasma 
phase, the free energy slowly approaches that of a free gluon gas at high temperature, 
and the speed of sound starts from a small value at T c and approaches the confor- 
mal value c 2 s = 1/3 as the temperature increases. The deviation from conformal 
invariance is strongest at T c , and is signaled by the presence of a non-trivial gluon 
condensate, which on the gravity side emerges as a deviation of the scalar solution 
that behaves asymptotically as r 4 close to the UV boundary. In the CP-violating 
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sector, the topological vacuum density TtFF has zero expectation value in the de- 
confined phase, in agreement with lattice results |21] and large- N c expectations. 

The analysis performed in |18| was completely general and did not rely on any 
specific form of the dilaton potential ^(A). In this paper we present instead a de- 
tailed analysis of an explicit model, whose thermodynamics matches quantitatively 
the thermodynamics of pure Yang-Mills theory. The (dimensionless) free energy, 
entropy density, latent heat and speed of sound, obtained on the gravity side by nu- 
merical integration of the 5D field equations, are compared with the corresponding 
quantities, calculated on the lattice for pure Yang-Mills at finite-T, resulting in ex- 
cellent agreement, for the temperature range that is accessible by lattice techniques. 
The same model also shows a good agreement with the lattice calculation of glueball 
mass ratios at zero temperature, and we find that the value of the deconfining critical 
temperature (in units of the lowest glueball mass) is also in good agreement with the 
lattice results. 

In short, the model we present gives a good phenomenological holographic de- 
scription of most static properties 1 (spectrum and equilibrium thermodynamics) of 
large- N c pure Yang-Mills, as computed on the lattice, for energies up to several times 
T c . Thus it constitutes a good starting point for the computation of dynamical ob- 
servables in a realistic holographic dual to QCD (as opposed to e.g. Af = 4 SYM), 
such as transport coefficients and other hydrodynamic properties that are not easily 
accessible by lattice techniques, at energies and temperatures relevant for relativistic 
heavy-ion collision experiments. We will report on such a calculation in the near 
future. 

The vacuum solution in this model is described in terms of two basic bulk fields, 
the metric and the dilaton. These are not the only bulk fields however, as the bulk 
theory is expected to have an a priori infinite number of fields, dual to all possible 
YM operators. In particular we know from the string theory side that there are a few 
other low mass fields, namely the RR axion (dual to the QCD #-angle) the NSNS and 
RR two forms B 2 and C 2 as well as other higher-level fields. With the exception of 
the RR axion, such fields are dual to higher- dimension and/or higher-spin operators 
of YM. Again, with the exception of the RR axion, they are not expected to play 
an important role into the structure of the vacuum and this is why we neglect them 
when we solve the equations of motion. However, they are going to generate several 
new towers of glueball states beyond those that we discuss in this paper (namely the 
++ glueballs associated to dilaton fluctuations, 2 ++ glueballs associated to graviton 
fluctuations and 0~ + glueballs associated to RR axion fluctuations). Such fields can 
be included in the effective action and the associated glueball spectra calculated. 
Since we do not know the detailed structure of the associated string theory, their 
effective action will depend on more semi-phenomenological functions like Z(X) in 

1 There are very few observables also that are not in agreement with YM. They are discussed in 
detail in ill. 
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( |5.8|) . These functions can again be determined in a way similar to Z(X). In partic- 
ular including the B 2 and C 2 field will provide l + ~ glueballs among others. Fields 
with spin greater than 2 are necessarily stringy in origin. We will not deal further 
with extra fields, like B 2 and C 2 and other as they are not particularly relevant for 
the purposes of this model, namely the study of finite temperature physics in the 
deconfined case. We will only consider the axion, as its physics is related to the 
CP-odd sector of YM with an obvious phenomenological importance. 

It is well documented that string theory duals of YM must have strong curvatures 
in the UV regime. This has been explained in detail in |15j where it was also ar- 
gued, that although the asymptotic AdS boundary geometry is due to the curvature 
non-linearities of the associated string theory, the inwards geometry is perturbative 
around AdS, with logarithmic corrections, generating the YM perturbation theory. 
The present model is constructed so that it takes the asymptotic AdS geometry for 
granted, by introducing the associated vacuum energy by hand, and simulates the 
perturbative YM expansion by an appropriate dilaton potential. In the IR, we do not 
expect strong curvatures in the string frame, and indeed the preferred backgrounds 
have this property. In this sense the model contains in itself the relevant expected 
effects that should arise from strong curvatures in all regimes. These issues have 
been explained in |14| and in more detail in f!5 |. 

The paper is organized a follows. In Section |2] we review the general features 
of the 5D model and its thermodynamics. In Section |] we discuss the issue of 
the possible sources of scheme dependence in the relation between the 5D dilaton 
field and the 4D Yang-Mills coupling. In Section [| we fix the form of the dilaton 
potential and we discuss in detail the independent parameters appearing in the model 
and the role they play in the dynamics. Section [| contains the main results: the 
comparison between our model and the lattice data for the thermodynamics of the 
deconfined phase, and for the spectrum of the confined phase. Finally, In Section 
H we provide some general conclusions and discuss future directions. Appendix [A| 
contains a description of the strategy we employ to numerically integrate Einstein's 
equations in order to find black hole solutions with different temperatures but same 
UV asymptotics. 



2. The 5D model 



The holographic dual of large N c Yang Mills theory, proposed in |f[4 |, is based on a 
five-dimensional Einstein-dilaton model, with the action: 



S 5 = -MlN 2 c I d°x^ 



i?--(<9$) 2 + U($) 



+ 2M p 3 iV 6 2 / d 4 xVhK. (2.1) 

JdM 



Here, M p is the five-dimensional Planck scale and N c is the number of colors. The 
last term is the Gibbons-Hawking term, with K being the extrinsic curvature of the 
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boundary. The effective five- dimensional Newton constant is G5 = l/(16nMpN 2 ), 
and it is small in the large- N c limit. 

Of the 5D coordinates {xj, r} i=0 ^ 3 , Xi are identified with the 4D space-time 
coordinates, whereas the radial coordinate r corresponds to the 4D RG scale. We 
identify A = with the running 't Hooft coupling X t = N c gy M , up to an a priori 
unknown multiplicative factor 2 , A = nX t . 

The dynamics is encoded in the dilaton potential 3 , V(X). The small-A and large- 
A asymptotics of V(X) determine the solution in the UV and the IR of the geometry 
respectively. For a detailed but concise description of the UV and IR properties of 



the solutions the reader is referred to Section 2 of [fil|. Here we will only mention 
the most relevant information: 

1. For small A, V^(A) is required to have a power-law expansion of the form: 

V(X) ~ ^(l + v X + Vl X 2 + ...), A^O. (2.2) 

The value at A = is constrained to be finite and positive, and sets the UV 
AdS scale t. The coefficients of the other terms in the expansion fix the /3- 
function coefficients for the running coupling X(E). If we identify the energy 
scale with the metric scale factor in the Einstein frame, as in |I4]| , we have: 

/3(A) = -^- = -M 2 -M 3 + ... 
a log h 

On = -Vn, Oi = —V\ V n . (2.3) 

° 8 ' 4 256 K J 



2. For large A, confinement and the absence of bad singularities 4 require: 
V(A)~A*(logA)' A^oo, {f^ 2 ^' Parbltra ^ 



(2.4) 



In particular, the values Q = 2/3, P =1/2 reproduce a linear glueball spec- 
trum, rn? n ~ n, besides confinement. We will restrict ourselves to this case in 
what follows. 

In the large N c limit, the canonical ensemble partition function of the model 
just described, can be approximated by a sum over saddle points, each given by a 
classical solution of the Einstein-dilaton field equations: 

Z{p) ~ e~ SlW + e~ S2W + ... (2.5) 

2 This relation is well motivated in the UV, although it may be modified at strong coupling (see 
Section |^) . The quantities we will calculate do not depend on the explicit relation between A and 
A*. 

3 With a slight abuse of notation we will denote V(X) the function V($>) expressed as a function 
of A = e*. 

4 We call "bad singularities" those that do not have a well defined spectral problem for the 
fluctuations without imposing extra boundary conditions. 
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where S{ are the euclidean actions evaluated on each classical solution with a fixed 
temperature T = 1//3, i.e. with euclidean time compactified on a circle of length 
0. There are two possible types of Euclidean solutions which preserve 3-dimensional 
rotational invariance. In conformal coordinates these are: 



1. Thermal gas solution, 

ds 2 = b 2 D (r) (dr 2 + dt 2 + dx m dx m ) , $ = $ (r), 
with r G (0, oo) for the values of P and Q we are using; 

2. Black hole solutions, 

dr 2 



(2.6) 



ds 2 = b(rY 



lf(r) 



f(r)dt 2 



doc Yfidx 



$ = $(r), 



(2.7) 



with r G (0,r h ), such that /(0) = 1, and /(77J = 0. 



In both cases Euclidean time is periodic with period (3 a and respectively for the 
thermal gas and black-hole solution, and 3-space is taken to be a torus with volume 
V 3o and V3 respectively, so that the black hole mass and entropy are finite 5 . 

The black holes are dual to a deconfined phase, since the string tension vanishes 
at the horizon, and the Polyakov loop has non- vanishing expectation value ( [p2[ p3|). 
On the other hand, the thermal gas background is confining. 

The thermodynamics of the deconfined phase is dual to the 5D black hole ther- 
modynamics. The free energy, defined as 



F=E-TS, 



(2.1 



is identified with the black hole on-shell action; as usual, the energy E and entropy S 
are identified with the black hole mass, and one fourth of the horizon area in Planck 
units, respectively. 

The thermal gas and black hole solutions with the same temperature differ at 
0(r 4 ): 



b(r) = b (r) 



1+ £77 + 



f(r) 



1 



Cr 4 



+ 



r -)■ 0, 



(2.9) 



where Q and C are constants with units of energy. As shown in [|18) they are related 
to enthalpy TS and the gluon condensate (TrF 2 ) : 



C 



TS 



22 (Tr F 2 ) T - (Tr F 2 ) a 



3(4vr) 



240M p 3 iV c 2 



(2.10) 



5 The periods and 3-space volumes of the thermal gas solution are related to the black-hole 
solution values by requiring that the geometry of the two solutions are the same on the (regulated) 
boundary. See p8| for details. 
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Although they appear as coefficients in the UV expansion, C and Q are determined 
by regularity at the black hole horizon. For T and S the relation is the usual one, 



T = -^M S=^- = 4n(M*N?V,)b*(r h ). (2.11) 

For Q the relation with the horizon quantities is more complicated and cannot be 
put in a simple analytic form. However, as discussed in p8| , for each temperature 



there exist only specific values of Q (each corresponding to a different black hole) 
such that the horizon is regular. 

At any given temperature there can be one or more solutions: the thermal gas 
is always present, and there can be different black holes with the same temperature. 
The solution that dominates the partition function at a certain T is the one with 
smallest free energy. The free energy difference between the black hole and thermal 



gas was calculated in |L8[ to be: 



150 -T- (2-12) 



M*N*V 3 M^N^V 3 " 4 

For a dilaton potential corresponding to a confining theory, like the one we will 
assume, the phase structure is the following 

1. There exists a minimum temperature T min below which the only solution is the 
thermal gas. 

2. Two branches of black holes ("big" and "small") appear for T > T min , but 
the ensemble is still dominated by the confined phase up to a temperature 
T > T ■ 

- L c " ± mm 

3. At T = T c there is a first order phase transition to the black hole phase. The 
system remains in the black hole (deconfined) phase for all T > T c . 

In principle there could be more than two black hole branches, but this will not 
happen with the specific potential we will use. 



3. Scheme dependence 

There are several sources of scheme dependence in any attempt to solve a QFT. Dif- 
ferent parametrizations of the coupling constant (here A) give different descriptions. 
However, physical statements must be invariant under such a change. In our case, 
reparametrizations of the coupling constant are equivalent to radial diffeomorphisms 
as we could use A as the radial coordinate. 

In the holographic context, scheme dependence related to coupling redefinitions 
translates into field redefinitions for the bulk fields. As the bulk theory is on-shell, 
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all on-shell observables (that are evaluated at the single boundary of space-time) are 
independent of the field redefinitions showing that scheme-independence is expected. 
Invariance under radial reparametrizations of scalar bulk invariants is equivalent 
to RG invariance. Because of renormalization effects, the boundary is typically 
shifted and in this case field redefinitions must be combined with appropriate radial 
diffeomorphisms that amount to RG-transformations. 

Another source of scheme dependence in our setup comes from the choice of 
the energy function. Again we may also consider this as a radial coordinate and 
therefore it is subject to coordinate transformations. A relation between A and E is 
the j3 function, 

/3(A). (3.1) 



dlogE 

(3 by definition transforms as a vector under A reparametrizations and as a form 
under E reparametrizations. /3(A) can therefore be thought of as a vector field 
implementing the change of coordinates from A to E and vice-versa. 

Physical quantities should be independent of scheme. They are quantities that 
are fully diffeomorphism invariant. If the gravitational theory had no boundary 
there would be no diffeomorphism invariant quantities, except for possible topological 
invariants. Since we have a boundary, diffeomorphism invariant quantities are defined 
at the boundary. 

Note that scalar quantities are not invariant. To be invariant they must be scalar 
and constant. We therefore need to construct scalar functions that are invariant 
under changes or radial coordinates. 

We can fix this reparametrization invariance by picking a very special frame. For 
example choosing the (string) metric in the conformal frame 

ds 2 = e 2A [dr 2 + dx^dx^ , A(r) (3.2) 

or in the domain-wall frame 

ds 2 = du 2 + e 2A dx^dx^ , X(u) (3.3) 

fixes the radial reparametrizations almost completely. In conformal frame, common 
scalings of r, x^ are allowed, corresponding to constant shifts of A(r). 

Eventually we are led to calculate and compare our results to other ways of 
calculating (like the lattice). Some outputs are easier to compare (for example cor- 
relators). Others are much harder as they are not invariant (like the value of the 
coupling at a given energy scale). 

In the UV such questions are well understood. The asymptotic energy scale is 
fixed by comparison to conformal field theory examples. This is possible because the 
space is asymptotically AdSs 6 . 



3 As the dilaton is now not constant there is a non-trivial question: in which frame is the metric 
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The coupling constant is also fixed to leading order from the coupling of the 
dilaton to D3 branes (up to an overall multiplicative factor). Subleading (in pertur- 
bation theory) redefinitions of the coupling constant and the energy lead to changes 
in the /3-function beyond two loops. 

More in detail, as it has been described in [[14], [OJ, the general form of the kinetic 
term for the gauge fields on a D 3 brane is expected to be: 

S F 2=e-*Z(R,0Tr[F 2 } , £= -e 2 ^ (3.4) 

where Z(R,£) is an (unknown) function of curvature R and the five-form field 
strength, £. At weak background fields, Z ~ — | + • • • . In the UV regime, expanding 



near the boundary in powers of the coupling A = N c e <t> we obtain, [15 



S F2 = N c Tr[F 2 ] ~ 



(3.5) 



where F(R,£) is the bulk effective action and i?*, are the boundary values for these 
parameters. Therefore the true 't Hooft coupling of QCD is 



A't Hooft 



Z (i?*, £*) 



Z ( R* , ) ( R* , ) ^ 



(3.6) 



In the IR, more important changes can appear between our A and other definitions 
as for example in lattice calculations. 

In the region of strong coupling we know much less in order to be guided con- 
cerning the correct definition of the energy. We can obtain some hints however by 
comparing with lattice results. 7 In particular, based on lattice calculations using the 
Schrondiger functional approach , it is argued that at long distance L the 't Hooft 
coupling constant scales as 

A lat ~ e mL , m ~ ^m 0++ . (3.7) 

This was based on a specific definition of the coupling constant, and length scale on 
the lattice as well as on numerical data, and some general expectations on the fall-off 
of correlations in a massive theory. This suggests an IR function of the form 

^ = AlogA , A = Aq e mL . (3.8) 



AdS. In [JFJ] it was argued that this should be the case in the string frame. The difference of course 
between the string and Einstein frame is subleading in the UV as the coupling constant vanishes 
logarithmically. But this may not be the case in the IR where we have very few criteria to check. 
In the model we are using we impose that the space is asymptotically AdS in the Einstein frame as 
this is the only choice consistent with the whole framework. 

7 We would like to thank K. Kajantie for asking the question, suggesting to compare with lattice 
data, and providing the appropriate references. 
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On the other hand our /3 function at strong coupling uses the UV definition of 
energy, log-E = Ae (the scale factor in the Einstein frame), E ~ 1/L and is 



dX _ 3 
L dL~ 2 A 



3a - 1 1 

1 + 



, A ~ 




(i) 







(3.9) 



4 a log A 

Consider now taking as length scale the string scale factor e As in the IR. 8 Since it 
increases, it is consistent to consider it as a monotonic function of length. From its 
relation to the Einstein scale factor A s = Ae + § log A and (|3.9|) we obtain 

dA 2a 



dA s a — 1 



AlogA+--- . (3.10) 



Therefore if we define as length scale in the IR 

2a a 2a 

\ogL = - — -A s L=(e As ) a - 1 (3.11) 

we obtain a running of the coupling compatible with the given lattice scheme. Note 
however that L = (e ) a_1 cannot be a global choice but should be only valid in the 
IR. The reason is that this function is not globally monotonic. 

We conclude this section by restating that physical observables are independent 
of scheme. But observables like the 't Hooft coupling constant do depend on schemes, 
and it is obvious that our scheme is very different from lattice schemes in the IR. 



4. The potential and the parameters of the model 

We will make the following ansatz for the potential, 

V(X) = p {l + V X + V^ 3 [log (1 + V 2 X 4/3 + V 3 A 2 )] V2 } , (4.1) 

which interpolates between the two asymptotic behaviors ( [2.2|) for small A and (|2.4j ) 
for large A, with Q = 2/3 and P = 1/2. Not all the parameters entering this potential 
have physical relevance. Below we will discuss the independent parameters of the 
model, and their physical meaning. 

The normalization of the coupling constant A. As discussed in the previous 
section, the relation between the bulk field A(r) and the physical QCD 't Hooft 
coupling X t = g\ M N c is a priori unknown. In the UV, the identification of the D3- 
brane coupling to the dilaton implies that the relation is linear, and depends on an 
a priori unknown coefficient k, defined as: 

A = KX t . (4.2) 

8 The string scale factor is not a monotonic function on the whole manifold, |l4| and this is the 
reason that it was not taken as a global energy scale. In particular in the UV, e As decreases until it 
reaches a minimum. The existence of the minimum is crucial for confinement. After this minimum 
e As increases and diverges at the IR singularity. 
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The coefficient k can in principle be identified by relating the perturbative UV ex- 
pansion of the Yang-Mills /3-function, to the holographic /3-function for the bulk field 
A: 

22 ^1 

0(X t ) = -AA? - X X\ +... (3 = 3—^, ft = — Pi . . . (4.3) 

Q Q 207 

/3(A) = -6 A 2 - M 3 + ..., b = ft*, h = - Vl - ^-vl .... (4.4) 

8 4 zoo 

The two expressions ( [4.3| ) and ( |4.4| ) are consistent with a linear relation as in 

(f4.2p, and expanding the identity K(3 t (X t ) = (3(K\ t ) to lowest order leads to: 

« = Po/bo. (4.5) 

Therefore, to relate the bulk field A to the true coupling At one looks at the linear term 
in the expansion of the potential. More generally, the other /3 function coefficients are 
related by f3 n = K n+1 b n , and the combinations = (3 n / (3q +1 are ^-independent 

(however they are scheme- dependent for n > 2). 

As discussed in Section ^ the introduction of the coefficient k amounts to a 
field redefinition and therefore its precise value does not affect physical (scheme- 
independent) quantities. In this sense, k is not a parameter that can be fixed by 
matching some observable computed in the theory. Assuming the validity of the 
relation Q4.2|) , we could eventually fix k by matching a RG-invariant (but scheme- 
dependent) quantity, e.g. A at a given energy scale. 

However, as we discuss later in this section, rescaling A in the potential (thus 
changing k) affects other parameters in the models, that are defined in the string 
frame, e.g. the fundamental string length £ s : if we hold the physical QCD string 
tension fixed, the ratio (£ s /£) scales with degree —2/3 under a rescaling of k. 

An important point to keep in mind, is that the simple linear relation ( |4.2j ) 
may be modified at strong coupling, but again this does not have any effect on 
physical observables. As long as we compute RG-invariant and scheme-independent 
quantities, knowledge of the exact relationship X = F(X t ) is unnecessary. 

The AdS scale £. This is set by the overall normalization of the potential, and 
its choice is equivalent to fixing the unit of energy. It does not enter dimensionless 
physical quantities. As usual the AdS length at large N c is much larger than the 
Planck length (£ p ~ l/(M p Nc^ 3 ), independently of the 't Hooft coupling. 

The UV expansion coefficients of V(X). They can be fixed order by order by 
matching the Yang-Mills /3-function. We impose this matching up to two-loops in 
the perturbative expansion, i.e. 0(A 3 ) in /3(A). One could go to higher orders by 
adding additional powers of A inside the logarithm, but since our purpose is not to 
give an accurate description of the theory in the UV, we choose not to introduce 
extra parameters 9 . 

9 Moreover, higher order /3-function coefficients are known to be scheme-dependent. 



- 11 - 



Identifying the energy scale with the Einstein frame scale factor, log E = log b(r), 
we have the relation ( |4.4|) between the /3-function coefficients and the expansion 
parameters of V(X), with 

v = V , v l = V l ^/V 2 . (4.6) 

The term proportional to V2 in eq. ( |4.1| ) is needed to reproduce the correct value of 
the quantity bi/b^ = /3i//3q = 51/121, which is invariant under rescaling of A. Thus, 
V2 is not a free parameter, but is fixed in terms of Vq and V\ by: 



4 / 23 + 36^/g V 9 ^ _ 51 

^"^l 81^ ) ' 6 °-8 y °' bj- 121" (4 ' 7) 

As explained earlier in this section, when discussing the normalization of the 
coupling, fixing the coefficient Vq is the same as fixing the normalization k through eq. 
( |4.5|) . As we argued, the actual value of k should not have any physical consequences, 
so it is tempting to set V = 1 by a field redefinition, A — > X/V and eliminate this 
parameter altogether. 

In fact, most of the quantities we will compute are not sensitive to the value of 
V , but for certain quantities, such as the string tension, some extra care is needed. 
In general, we can ask whether two models of the same form ( |2.1| ), but with different 
potentials V(X) and V^(A), such that V^(A) = V(a\) for some constant a, lead to 
different physical predictions. As we can change from one model to the other simply 
by a field redefinition A — > aX ( this has no effect on the other terms in the action 
in the Einstein frame, eq. (|2.1| ) ), clearly the two potentials lead to the same result 
for any physical quantity that can be computed unambiguously from the Einstein 
frame action, e.g. dimensionless ratios between glueball masses, critical temperature, 
latent heat etc. 

However a rescaling of A does affect the string frame metric, since the latter 
explicitly contains factors of A: b s (r) = 5(r)A 2 / 3 thus, under the rescaling A — > 
a A, b s (r) — > a 2 ^ 3 b s (r). This means that any dimensionless ratio of two quantities, 
such that one of them remains fixed in the string frame and the other in the Einstein 
frame, will depend on a. An example of this is the ratio £ s /£, where £ s is the string 
length, that we will discuss shortly. 

Therefore, we can safely perform a field redefinition and set Vq to a given value, 
as long as we are careful when computing quantities that depend explicitly on the 
fundamental string length. 

Bearing this caveat in mind, we will choose a normalization such that bo = (3q, 

i.e. 

Vq = -A>, (4-8) 

so that the normalization of A in the UV matches the physical Yang-Mills coupling. 
With this choice, out of the four free parameters Vi appearing in ( |4.1| ) only V\ and 
V3 play a non-trivial role (V 2 being fixed by eq. ([4. 7|) ) . 



- 12 - 



The 5D Planck scale M p . M p appears in the overall normalization of the 5D 
action (|2.1|) . Therefore it enters the overall scale of quantities derived by evaluating 
the on-shell action, e.g. the free energy and the black hole mass. It also sets the 
conversion factor between the entropy and horizon area. M p cannot be fixed directly 
as we lack a detailed underlining string theory for YM. To obtain quantitative pre- 
dictions, M p must be fixed in terms of the other dimension-full quantity of the model, 
namely the AdS scale £. As shown in |T§j this can be done by imposing that the 



high-temperature limit of the black hole free energy be that of a free gluon gas with 
the correct number of degrees of freedom 10 . This requires: 

(M p £) 3 = JL. (4.9) 

The string length. In the non-critical approach the relation between the string 
length £ s and the 5D Planck length (or the AdS length £) is not known from first 
principles. The string length does not appear explicitly in the 2-derivative action 
( |2.1|) , but it enters quantities like the static quark- ant iquark potential. The ratio 
£ s /£ can be fixed phenomenologically to match the lattice results for the confining 
string tension. 

More in detail, the relation between the fundamental and the confining string 
tensions Tf and a is given by: 



a 



T f tf(r,)\*' 3 (r.), (4.10) 



where r* is the point where the string frame scale factor, b s (r) = 6(r)A 2//3 (r), has its 
minimum. Fixing the confining string tension by comparison with the lattice result 
we can find Tf (more precisely, the dimensionless quantity Tf£ 2 , since the overall scale 



of the metric depends on £). The string length is in turn given by £ s /£ = 1/ ■J 2irTf£ 2 . 
As is clear from eq. (|4.10 ), rescaling A — > a A, keeping the value of the QCD 



string tension a and of the AdS scale £ fixed, affects the fundamental string length 
in AdS units as £ s /£ — > a~ 2 / 3 (£ s /£). Therefore two models a and b, defined in the 
Einstein frame by eq. (|Q|), but with potentials related by Vb(X) = V a (a\), must 
have different fundamental string tensions in order to reproduce the same result for 
the QCD string tension. The quantity £ s /£ therefore depends on the value of V . 



Integration constants. Besides the parameters appearing directly in the gravita- 
tional action, there are also other physically relevant quantities that label different 
solutions to the 5-th order system of field equations. Any solution is characterized by 

10 Note that this is conceptually different from the J\f = 4 case. There, near the boundary, the 
theory is strongly coupled and this number must be calculated in string theory. It is different by a 
factor of 3/4 from the free sYM answer. Here near the boundary the theory is free. Therefore the 
number of degrees of freedom can be directly inferred. 
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a scale A, the temperature T and a value for the gluon condensate Q, that correspond 
to three of the five independent integration constants. 11 

Regularity at the horizon fixes Q as a function of T, so that effectively the gluon 
condensate is a temperature-dependent quantity. 

The quantity A controls the asymptotic form of the solution, as it enters the dila- 
ton running in the UV: A ~ -(b Q logrA)" 1 . It can be defined in a reparametrization 
invariant way as: 



cxp 



j_ 

b \ 



A = t~ x lim < 6(A) h^}, (4-11) 

and it is fixed once we specify the value of the scale factor b(X) at a given Ao- 

Every choice of A corresponds to an inequivalent class of solutions, that differ by 
UV boundary conditions. Each class is thermodynamically isolated, since solutions 
with different A's have infinite action difference. Thus, in the canonical partition sum 
we need to consider only solutions with a fixed value of A. However, this choice is 
merely a choice of scale, as solutions with different A's will give the same predictions 
for any dimensionless quantity. In short, A is the holographic dual to the QCD strong 
coupling scale: it is defined by the initial condition to the holographic RG equations, 
and does not affect dimensionless quantities such as mass ratios, etc. Therefore, as 
long as all solutions we consider obey the same UV asymptotics, the actual value 
of A is immaterial, since the physical units of the system can always be set by fixing £. 

To summarize, the only nontrivial phenomenological parameters we have at our 
disposal are V\ and V3 appearing in (|4.1|). The other quantities that enter our model 



are either fixed by the arguments presented in this section, or they only affect trivially 
(e.g. by overall rescaling that can be absorbed in the definition of the fundamental 
string scale) the physical quantities. 

In the next section we present a numerical analysis of the solutions and thermo- 
dynamics of the model defined by eq. (|4.1|) , and show that for an appropriate choice 
of the parameters it reproduces the lattice results for the Yang-Mills deconfinement 
transition and high-temperature phase as well as the zero temperature glueball data. 

5. Matching the thermodynamics of large- N c YM 



Assuming a potential of the form ([O]), we look for values of the parameters such 
that the thermodynamics of the 5D model match the lattice results for the thermo- 



dynamics of 4D YM. As explained in Section ||, we set V and V 2 as in eqs. ([18]) and 
( ET7D , respectively, with b = /3 = 22/3 (4vr)- 2 . 



11 The remaining two are the value /(0) which should be set to one for the solution (^7) to obey 
the right UV asymptotics, and an unphysical degree of freedom in the reparametrization of the 
radial coordinate. 
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We then vary V\ and V3 only. We fix these parameters by looking at thermody- 
namic quantities corresponding to the latent heat per unit volume, and the pressure 
at one value of the temperature above the transition, which we take as 2T C . 

It is worth remarking that Vy and V3 are phenomenological parameters that we 
use to fit dimensionless QCD quantities. The single (dimension-full) parameter of 
pure Yang-Mills, the strong coupling scale, is an extra input that fixes the overall 
energy scale of our solution. 

Using the numerical method explained in the Appendix, for each set of parame- 
ters (Vi, V3) we numerically generate black hole solutions for a range of values of A^, 
then from the metric at the horizon and its derivative we extract the temperature 
and entropy functions T(X h ) and S(Xh), and the function ^(A^) from the integrated 
form of the first law, eq (|A.10|) . The behavior of these functions is shown in Figures 
H, H and f|, for the best fit parameter values that we discuss below. One can see the 



existence of a minimal temperature T„ 



and a critical value A r where J 7 



changes sign. The resulting function JF(T) is shown in Figure [|. 




-0.02 



-0.03 



Figure 1: The Free energy density (in units of T c ) as a function of T/T c , for V\ = 14 
and V3 = 170. The vertical lines correspond to the critical temperature (solid) and the 
minimum black hole temperature (dashed). 

The phase transition is first order, and the latent heat per unit volume L^, 
normalized by N^T^, is given by the derivative of the curve in Fig. [I] at T/T c = 1. 
Equivalently, Lh is proportional to the jump in the entropy density s = S/V3 at 
the phase transition from the thermal gas (whose entropy is of 0(1), in the limit 
iV c — > 00) to the black hole (whose entropy scales as N% in the same limit): thus, in 
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the large N c limit, 

Lh = TAs ~ T c s(A c ) (5.1) 

up to terms of 0(l/iV 2 ). 

At this point we can look for values of V\ and V 3 that best fit the lattice data for 
the deconfined phase of thermal Yang-Mills. We compare our results to the data of 
G. Boyd et al. ||25|| . The relevant quantities to compare are the dimensionless ratios 
p(T)/T 4 , e(T)/T 4 and s(T)/T 3 , where p = T/V3 is the pressure, and e = p + Ts 
is the energy density. Lattice results for these functions are available in the range 
T = T c ~ 5T C , and can be seen in Figure 7 of f25| . The analysis of [^5| correspond to 



N c = 3, but one expects that the thermodynamic functions do not change to much 
for large N c 12 . 

An additional quantity of relevance is the value for the "dimensionless" latent 
heat per unit volume, Lh/T^ which for large N c was found in |27|] to be (Lh/T^iat = 
0.31iV 2 . The result for iV 2 = 3 is slightly lower ( ~ 0.28A" 2 ). 

As already noted in |L7|, the qualitative features of the thermodynamic func- 
tions are generically reproduced in our setup: the curves 3p(T)/T 4 , e(T)/T 4 and 
3s(T)/4T 3 increase starting at T c , then (very slowly) approach the constant free 
field value 7r 2 iV 2 /15 (given by the Stefan-Boltzmann law) as T increases. By com- 
puting the thermodynamic functions for various sets of values of V\ and V3 we find 
that 1) V\ roughly controls the height reached by the curves p(T)/T 4 , e(T)/T 4 and 
s(T)/T 3 at large T/T c (~ a few): for larger V\ the curves approach the free field 
limit faster; 2) V3 does not affect much the height of the curves at large T/T c , but 
on the other hand it changes the latent heat, which is increasing as V3 decreases. 

We find that the best fit corresponds to the values 

Vi = 14 V 3 = 170. (5.2) 

Below we discuss the values of various physical quantities ( both related to thermody- 
namics, and to zero-temperature properties) obtained with this choice of parameters. 

5.1 Latent heat and equation of state 

The comparison between the curves p(T)/T 4 , e(T)/T 4 and s(T)/T 3 obtained in our 
models with ( p.2[ ), and the lattice results [^SJ 13 is shown in Figure [5[ The match is 
remarkably good for T c < T < 2T C , and deviates slightly from the lattice data in the 
range up to 5T C . 

The latent heat we obtain is: 

L ft /T 4 = 0.31iV 2 , (5.3) 
which matches the lattice result for N c — > 00 \2%. 



12 See e.g. |2(J, in which results for iV c = 8 do not different significantly from those for N c 
We thank B. Lucini for useful correspondence on this point. 
13 We thank F. Karsch for providing us the relevant data. 
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Figure 2: Temperature in units of T m i n , as a function of A^, for V\ = 14 and V3 = 170. 
The dashed horizontal and vertical lines indicate the critical temperature and the critical 
value of the dilaton field at at the horizon 
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Figure 3: The free energy density in units of T m j n , as a function of A^ 

An interesting quantity is the trace anomaly (e — 3p)/T 4 , (also known as interac- 
tion measure), that indicates the deviation from conformality, and it is proportional 
to the gluon condensate. The trace anomaly in our setup is shown, together with 
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Figure 4: Entropy density in units of T m j n , as a function of 

the corresponding lattice data, in figure (f|), and the agreement is again very good. 

We also compute the specific heat per unit volume c v , and the speed of sound c s 
in the deconfined phase, by the relations 

r^d 2 F , s , r t . 

These are shown in Figures |8| and [7] respectively. The speed of sound is shown 
together with the lattice data, and the agreement is remarkable. 

5.2 Glueball spectrum 



In JL4J], the single phenomenological parameters of the potential was fixed by looking 
at the zero-temperature spectrum, i.e. by computing various glueball mass ratios and 
comparing them to the corresponding lattice results. The associated thermodynamics 



for this potential was studied in Jl7| which was in qualitative agreement with lattice 
QCD results, but not in full quantitative agreement. This is due to the fact that 
the thermodynamics depends more on the details of the potential than the glueball 
spectrum for the main Regge trajectories. Here we use the potential (|4.1|), but 
with the two phenomenological parameters V\ and V3 already determined by the 
thermodynamics (|5.2| ). 

The glueball spectrum is obtained holographically as the spectrum of normal- 
izable fluctuations around the zero-temperature background. As explained in the 
introduction, and motivated in [FJ, |T3|, here we consider explicitly the 5D metric, 
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Figure 5: Temperature dependence of the dimensionless thermodynamic densities s/T 3 
(light blue), p/T 4 (dark blue) and e/T (green), normalized such that they reach the 
common limiting value 7r 2 /15 (dashed horizontal line) asT-f oo. The dots correspond to 
the lattice data for N c = 3 pq ]. 



one scalar field (the dilaton), and one pseudoscalar field (the axion). As a conse- 
quence, the only normalizable fluctuations above the vacuum correspond to spin 
and spin 2 glueballs 14 (more precisely, states with J PC = ++ , CT + , 2 ++ ), each species 
containing an infinite discrete tower of excited states. 

In 4D YM there are many more operators generating glueballs, corresponding 
to different values of J PC , that are not considered here. These are expected to cor- 
respond holographically to other fields in the noncritical string spectrum (e.g. form 
fields, which may yield spin 1 and CP-odd spin 2 states) and to higher string states 
that provide higher spin glueballs. As the main focus of this paper is reproducing the 
YM thermodynamics in detail rather than the entire glueball spectrum, we chose not 
to include these states 15 . Therefore we only compare the mass spectrum obtained 
in our model to the lattice results for the lowest ++ ,0 h ,2 ++ glueballs and their 
available excited states. These are limited to one for each spin species, and none 



for the spin 2, in the study of [p8| , which is the one we use for our comparison. This 



14 Spin 1 excitations of the metric can be shown to be non-normalizable. 

15 A hither reason is that, unlike the scalar and (to some extent) the pseudoscalar sector that we 
are considering, the action governing the higher Regge slopes is less and less universal as one goes 
to higher masses. Only a precise knowledge of the underline string theory is expected to provide 
detailed information for such states. 
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Figure 6: The trace anomaly as a function of temperature in the deconfined phase of the 
holographic model (solid line) and the corresponding lattice data [25| for N c = 3 (dots). 
The peak in the lattice data slightly above T c is expected to be an artifact of the finite 
lattice volume. In the infinite volume limit the maximum value of the curve is at T c , and 
it equals L h /N?T*. 





HQCD 


iV c = 3| 


jV c = oo§ 


m *++/m ++ 


1.61 


1.56(11) 


1.90(17) 


m 2 ++ / m ++ 


1.36 


1.40(4) 


1.46(11) 



Table 1: Glueball Masses 



provides two mass ratios in the CP-even sector and two in the CP-odd sector. 

The glueball masses are computed by first solving numerically the zero-temperature 
Einstein's equations, obtained from ( A.l ) by setting f(r) = 1, and using the resulting 
metric and dilaton to setup an analogous Schrodinger problem for the fluctuations, 
14]. The results for the parity-conserving sector are shown in Table [l], and are in 



good agreement with those reported by p8[ for N c = 3, whereas the results reported 
by f29| for large iV c are somewhat larger. The CP-violating sector (axial glueballs) 
will be discussed separately. 

We should add that there are other lattice studies (see e.g. PTJ) that report 
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Figure 7: The speed of sound in the deconfined phase, as a function of temperature, for 
the holographic model (solid line) and the corresponding lattice data [25 for N c = 3 (dots). 
The dashed horizontal line indicates the conformal limit (? s = 1/3. 



additional excited states. Our mass ratios offer a somewhat worse fit of the mass 
ratios found in |30[ (whose results are not entirely compatible with those of [28 



for the states the two studies have in common). We should stress however that 
reproducing the detailed glueball spectrum is secondary here since the main focus is 
thermodynamics. However, the comparison of our spectrum to the existing lattice 
results shows that our model provides a good global fit to 4D YM also with respect 
to quantities beyond thermodynamics. 

Unlike the various mass ratios, the value of any given mass in AdS-\ength units 
(e.g. m 0++ £) does depend on the choice of integration constants in the UV, i.e. on 
the value of bjjy and Xuv- Therefore its numerical value does not have an intrinsic 
meaning. However it can be used as a benchmark against which all other dimension- 
full quantities can be measured (provided one always uses the same UV b.c. ). On 
the other hand, given a fixed set of initial conditions, asking that mo++ matches the 
physical value (in MeV) obtained on the lattice, fixes the value of t hence the energy 
unit. 
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Figure 8: The specific heat (divided by T 3 ), as a function of temperature, in the deconfined 
phase of the holographic model. 

5.3 Critical Temperature 

The thermodynamic quantities we have discussed so far, are dimensionless ratios, 
in units of the critical temperature. To compute T c , we need an extra dimension- 
full quantity which can be used independently to set the unit of energy. In lattice 
studies this is typically the confining string tension a in the T = vacuum, with 
a value of around (440MeV) 2 , and results are given in terms of the dimensionless 
ratio T c /y/a. In our case we cannot compute a directly, since it depends on the 
fundamental string tension, which is a priori unknown. Instead, we take the mass 
m 0++ of the lowest-lying glueball state as a reference. 

with V\ and V3 fixed as in 



We compute m ++ with the potential (fO 



then compare T c /m 0++ to the same quantity obtained on the lattice. For the lattice 
result, we take the large N c result of [p7| , T c /yfa = 0.5970(38), and combine it with 
the large N c result for the lowest-lying glueball mass |p9fl , m 0++ /^ = 3.37(15). The 
two results are in fair agreement, without need to adjust any extra parameter: 

— J =0.167, (— ] =0.177(7). 

™QJhQCD \ m J lattice 

In physical units, the critical temperature we obtain is given by 



(5.5) 



T c = 0.56 v/a = 247 Me V. 



(5.6) 
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5.4 String tension 

The fundamental string tension Tf = —72 cannot be computed from first principles in 
our model, but can be obtained using as extra input the lattice value of the confining 
string tension a, at T = 0. The fundamental and confining string tensions are related 
by eq. ( 4.10|) . As for the critical temperature, we can relate Tf to the to the value 



of the lowest-lying glueball mass, by using the lattice relation y/a = m 0++ /3.37 



29]. Since what we actually compute numerically is m 0++ £, this allows us to obtain 



the string tension Tf (and fundamental string length l s = 1/ yj2nTf) in AdS units. 
Evaluating numerically the factor ^(r^A 2 / 3 ^*) that appears in eq. ( f4.10| ), we find: 



7>r = 6.5, e B /£ = 0.15. (5.7) 

This shows that the fundamental string length in our model, although not paramet- 
rically small, turns out to be one order of magnitude smaller than the AdS length, 
which is a good sign for the validity of the derivative expansion. The meaning of 
this fact is a little more complicated conceptually, as the discussion in |T5| indicates. 



Also, we should stress that, as discussed in Section [|, this result depends on our 
choice of the overall normalization of A: changing the potential by A — > kX will yield 
different numerical values in ( |5.7|) without affecting the other physical quantities. 

5.5 CP-odd sector 

The CP-odd sector of pure Yang-Mills is described holographically by the addition 
of a bulk pseudoscalar field a(r) (the axion) with action: 16 

M 3 r 

S axion = / d 5 xZ(X)^(d^a)(d fl a) . (5.8) 



The field a(r) is dual to the topological density operator TtFF. The prefactor Z(X) 
is a dilaton-dependent normalization. The axion action is suppressed by a factor 
1/iV 2 with respect to the action (|2 . 1| ) for the dilaton and the metric, meaning that in 
the large- N c limit one can neglect the back-reaction of the axion on the background. 



As shown in [II], requiring the correct scaling of a(r) in the UV, and phe- 
nomenologically consistent axial glueball masses, constrain the asymptotics of Z(X) 
as follows: 

Z{\) ~ Z Q , A ->■ 0; Z{\) ~ A 4 , A -> oo, (5.9) 

where Z is a constant. As a simple interpolating function between these large- and 
small- A asymptotics we can take the following: 

Z(X) = Z (l + c a X 4 ). (5.10) 



16 This action was justified in |l5| . The dilaton dependent coefficient Z(X) is encoding both 
the dilaton dependence as well as the UV curvature dependence of the axion kinetic terms in the 
associated string theory. We cannot determine it directly from the string theory, but we pin it down 
by a combination of first principles and lattice input, as we explain further below. 
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The parameter Zq can be fixed by matching the topological susceptibility of pure 
Yang-Mills theory, whereas c a can be fixed by looking at the axial glueball mass 
spectrum. 



Axial glueballs. As in [[14), we can fix c a by matching to the lattice results the 

mass ratio tuq |_ / m 0++ between the lowest-lying axial and scalar glueball states. This 

is independent of the overall coefficient Z in ( |5.10| ). The lattice value m (_/m ++ = 

1.49 PS is obtained for: 



c a = 0.26. (5.11) 

With this choice, the mass of the first excited axial glueball state is in good agreement 
with the corresponding lattice result PS|: 



^z±i) =2.10 (!^z±l) =2.12(10). (5.12) 

m 0++ / hQC D V m 0++ / i at ti ce 

Topological Susceptibility. In pure Yang-Mills, the topological x susceptibility 
is defined by: 

E{6)= l - X e\ (5.13) 

where E{9) is the vacuum energy density in presence of a ^-parameter. E{6) can be 
computed holographically by solving for the axion profile a(r) on a given background, 
and evaluating the action ( |5.8| ) on-shell. 

In the deconfined phase, the axion profile is trivial, implying a vanishing topo- 
logical susceptibility ||18|| . This is in agreement with large- N c arguments and lattice 



results [|lj. 

In the low-temperature phase, the axion acquires a non-trivial profile, 

a(r) = a uv F(r) S jT _*_ 5T . (5.14, 

This profile is shown, for the case at hand, in Figure || where the axion is normalized 
to its UV value. 

The topological susceptibility is given by [|T4[ : 

-i 



X = M'FiO)- 1 = Ml 



dr 



r 



(5.15) 



where Z[f) = Z(X(r)). Evaluating this expression numerically with Z(X) as in 
( |5.10p , and c a = 0.26 (to match the axial glueball spectrum ), we can determine 



the coefficient Zq by looking at the lattice result for \- F° r N c = 3, [31| obtained 
X = (l91MeV) 4 , which requires Z Q = 133. 

In Table ^ we present a summary of the various physical quantities discussed in 
this section, as obtained in our holographic model, and their comparison with the 
lattice results for large N c (when available) and for iV c = 3. The quantities shown 
in the upper half of the table are the ones that were used to fix the free parameters 
(reported in the last column) of the holographic model. 
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Figure 9: Axion profile in the radial direction. The x-axis is taken to be the energy scale, 
E(r) = Eob(r), where the unit Eq is fixed to match the lowest glueball mass. 



5.6 Coupling normalization 

Finally, we can relate the field A(r) to the running 't Hooft coupling. All other 
quantities we have discussed so far are scheme- independent and RG-invariant. This 
is not the case for the identification of the physical YM 't Hooft coupling, which is 
scheme dependent. 

In the black hole phase we can take A^ = A(r/j) as a measure of the temperature- 
dependent coupling. In figure we show X h as a function of the temperature in the 
range T c to 5T C . 

As a reference, we may take the result of [p5 , that found g 2 (5T c ) ~ 1.5 for N c = 3, 



which translates to A t (5T c ) ~ 5. On the other hand, if we make the assumption that 
the identification X = X t is valid at all scales (not only in the UV), we find in our 
model Xt(5T c ) ~ 0.04 (see Figure |T0|) , i.e. a factor of 100 smaller than the lattice 
result. 

This discrepancy is almost certainly due to the identification ( |4.2D being very 
different from lattice at strong coupling. 



6. Discussion and Outlook 

The construction presented in this paper offers a holographic description of large- A^ c 
Yang-Mills theory that is both realistic and calculable, and in quite good agreement 
with a large number of lattice results both at zero and finite temperature. 
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Figure 10: The coupling at the horizon as a function of temperature in the range T C -5T C 



It is a phenomenological model; as such it is not directly associated to an explicit 
string theory construction. In this respect it is in the same class as the models based 
on pure AdS backgrounds (with hard or soft walls) [0, §, [32|, |33|, or on IR deforma- 
tions of the AdS metric In comparison to them however, our approach has 
the advantage that the dynamics responsible for strong coupling phenomena (such 
as confinement and phase transitions) is made explicit in the bulk description, and 
it is tied to the fact that the coupling constant depends on the energy scale and 
becomes large in the IR. This makes the model consistent and calculable, once the 
5D effective action is specified: the dynamics can be entirely derived from the bulk 
Einstein's equation. The emergence of an IR mass scale and the finite temperature 
phase structure are built-in: they need not be imposed by hand and do not suffer 
from ambiguities related to IR boundary conditions (as in hard wall models) or from 
inconsistencies in the laws of thermodynamics (as in non-dynamical soft wall models 
based on a fixed dilaton profile || or on a fixed metric |S5[ ) . More specifically, in 
our approach it is guaranteed that the Bekenstein-Hawking temperature of the black 
hole matches the entropy computed as the derivative of the free energy with respect 
to temperature. 

With an appropriate choice of the potential, we provide a realistic and quan- 
titatively accurate description of essentially all the static properties (spectrum and 
equilibrium thermodynamics) of the dynamics of pure Yang-Mills. The main ingredi- 
ent responsible for the dynamics (the dilaton potential) is fixed through comparison 
with both perturbative QCD and lattice results. It is worth stressing that such 
a matching on the quantitative level was only possible because the class of holo- 
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HQCD 


lattice N c 


= 3 


lattice iV c — > oo 


Parameter 


[p/(iV 2 T 4 )] T=2Tc 


1.2 


1.2 






VI = 14 


L h /{N 2 X) 


0.31 


0.28 


25 




0.31 \27\ 


V3 = 170 


Ip/(N?T*)} t ^ +00 


7T 2 /45 


7T 2 /45 




7T 2 /45 


M p £= [457T 2 ]- 1 / 3 


m ++ 1 J~o 

\J / V 


3.37 


3.56 


28 




3.37 129 


Ul = 0.15 


rriQ-+ /m Q ++ 


1.49 


1.49 


28 


] 




c a = 0.26 


X 


(191Me1/) 4 


(191Me\/) 4 


m 


- 


Z = 133 


T c /m ++ 


0.167 




0.177(7) 




m *++ / m ++ 


1.61 


1.56( 


11) 




1.90(17) 




m 2 ++/ m o++ 


1.36 


1.40(4) 




1.46(11) 




m *-+ /m ++ 


2.10 


2.12(10) 







Table 2: Collected in this table is the complete set of physical quantities that we computed 
in our model and compared with data. The upper half of the table contains the quantities 
that we used as input (shown in boldface) for the holographic QCD model (HQCD). Each 
quantity can be roughly associated to one parameter of the model (last column). The lower 
half of the table contains our "postdictions" (i.e. quantities that we computed after all the 
parameters were fixed) and the comparison with the corresponding lattice results. The 
value we find for the critical temperature corresponds to T c = 247MeV. 

graphic models we discuss generically provides a qualitatively accurate description of 
the strong Yang- Mills dynamics. This is a highly non-trivial fact, that strongly indi- 
cates that a realistic holographic description of real-world QCD might be ultimately 
possible. 

Although the asymptotics of our potential is dictated by general principles, we 
based our choice of parameters by comparing with the lattice results for the ther- 
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modynamics. There are other physical parameters in the 5D description that do not 
appear in the potential: the 5D Planck scale, that was fixed by matching the free 
field thermodynamics in the limit T — > oo; the coefficients in the axion kinetic term, 
that were set by matching the axial glueball spectrum and the topological suscepti- 
bility (from the lattice). The quantities that we use as input in our fit, as well as the 
corresponding parameters in the 5D model, are shown in the upper half of Table 0. 

The fact that our potential has effectively two free parameters depends on our 
choice of the functional form. This functional form contains some degree of arbi- 
trariness, in that only the UV and IR asymptotics of V(A) are fixed by general 
considerations (matching the perturbative /3-function in the UV, and a discrete lin- 
ear glueball spectrum for the IR). Therefore the results presented in this paper offer 
more a description, rather than a prediction of the thermodynamics. 

Nevertheless, there are several quantities that we successfully "postdict" (i.e. 
they agree with the lattice results) once the potential is fixed: apart from the good 
agreement of the thermodynamic functions over the whole range of temperature 
explored by the lattice studies (see Figures ||, |6| and 0), they are the lowest glueball 
mass ratios and the value of the critical temperature. The comparison of these 
quantities with the lattice results is shown in the lower half of Table |2], and one can 
see that the agreement is overall very good. Moreover the model predicts the masses 
of the full towers of glueball states in the + ~, ++ , 2 ++ families. 

The fact that our model is consistent with a large number of lattice results is 
clearly not the end of the story: its added value, and one of the main reasons for its 
interest lies in its immediate applicability beyond equilibrium thermodynamics, i.e. 
in the dynamic regimes tested in heavy-ion collision experiments. This is a generic 
feature of the holographic approach, in which there are no obstructions (as opposed to 
the lattice) to perform real-time computations and to calculate hydrodynamics and 
transport coefficients. Our model provides a framework to compute these quantities 
in a case where the static properties agree with the real- word QCD at the quantitative 
level, and this will be reported in a future publication. 

In the approach followed in this work, we took the phenomenological two- derivative 
bulk action at face value and applied the standard rules of holography to relate it 
to a 4D theory. Once the 5D action is taken as an input, this procedure is well de- 
fined, and yields unique and well defined results. The smallness of the gravitational 
coupling in the large N c limit ensures that quantum corrections can be neglected, 
independently of the details of the underlying fundamental theory. Here we did not 
discuss the justification of this approach, and its possible limits of validity, in the 
broader context of non-critical string theory. These are still open issues, that were 
addressed in [[H and more recently in |[L5|| . 
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APPENDIX 



A. Numerical technique 

To compute the thermodynamics from the 5D holographic model we need to find 
different temperature black hole solutions of Einstein-Dilaton field equations with the 
potential ( |4.1|) . This cannot be done analytically, and we must resort to numerical 
integration. The main difficulty in reaching this goal is how to fix the integration 
constant, or equivalently boundary conditions so that the following requirements are 
met: 

i. all the solutions must obey the same UV asymptotics as r — > 0. In particular 
f(r) should go to one a the the boundary and the integration constant A 
discussed in the Section £| should be the same for all solutions. 

ii. the solutions must have regular horizons, rather than naked singularities. As 
shown in ||18|| , this imposes one constrain among the values of b(r), b(r), f(r), 
and A(r) evaluated at the "would be" horizon, i.e. where / = 0. 

Clearly the first requirement privileges a numerical integration starting from the UV, 
while the second one starting from the horizon. One elegant way out of this dilemma 
is to use a solution-generating technique based on a set of symmetry transformations 
enjoyed by the field equations: we starts numerical integration at the horizon im- 
posing regular initial conditions, obtaining thus a regular solution; then we act with 
the symmetry transformations in order to generate a new solution that satisfy given 
UV boundary conditions. The second step can be done analytically, and since the 
symmetry transformation does not spoil regularity of the horizon, the new solution 
can satisfy both the above requirements. 

The equations of motion can be written in the following form 

A(r) = \\{r)sJb{r)W(r) , b(r) = -b(r)W(r) , (A.l) 
f{r) = 3b{r)f{r)W(r) , W(r) = Ab{r)W{rf - j$W(r) - $$V(\(r)) , 

where W(r) is an auxiliary function that we introduce for convenience and reduces 



to the zero-temperature superpotential in the thermal-gas solution 

The five integration constants discussed in the previous Section are fixed in the 
numerical procedure by giving initial conditions at the horizon and by exploiting 
the symmetries of the system of equations ( |A.1| ) to rescale the solutions at different 
temperatures so that they describe the same theory, i.e. with the same UV behavior 
(or, in other words, with the same 17 A). More precisely, we set the five initial 



17 The solutions at different temperature will also be characterized by the same normalization for 
/ in the UV and the same location of the UV boundary, rjjy. 
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conditions for A(r), b(r), f(r), f(r), W(r) at a point r = rj close to the horizon 
using the expansion of these functions to the first order in = (r^ — r«) (the dots 
in the formulae stand for O ((r^ — rj) 2 ) terms): 

Kn) = b h + b h e h + ... , X(n) = X h + X h e h + ... , (A.2) 
f(n) = f h e h + ... , W(r t ) = W h + W h e h + ... , 

where we have introduced the notation bh = b(rh) and similarly for the other dotted 
quantities appearing in the equations above. 

The solution is uniquely determined by specifying the five quantities {A^, bh, 
Wh, fh, £h) at the position r h 18 , since the linear (for small e/J system ( |A.2| ) can be 
translated to an initial value problem at rj with full set of five initial data, {A(rj), 

6(n), Win), f{n), f( n )}. 



Regularity at the horizon constrains the initial data to obey the relation [18 



W h = hhV{Xh) . (A.3) 

The other first order coefficients bh, Wh, Xh are all determined by solving the equations 
of motion close to the horizon. 

Among the (free) parameters characterizing the initial conditions Xh,bh, fh,^h, 
we vary the physical parameter Xh alone to scan solutions with different horizon value 
of A(r), while bh, fh, £h are kept fixed (and can be chosen arbitrarily). For each A^ we 
generate numerically a solution b(r), A(r), f(r). Apart for the value of A(r^), these 
these solutions all share the same near-horizon behavior, but for the time being they 
all have different UV asymptotics (including different values of /(0)). 

The second step of the numerical procedure is to match the solutions in the UV. 
In fact, in order for the solutions obtained for different A^'s to describe the same 
theory, we have to make sure that we have the same asymptotic behavior for the 
scale factor and the dilaton (i.e. the same A), and the same normalization for /. 
For convenience (but this is not a strict physical requirement), we can also ask that 
all the solutions have coinciding UV boundary, rather than coinciding horizon. For 
each solution obtained numerically, we can generate a new solution obeying these 
three requirements by acting with the following three transformations (parametrized 
by the constants Sf, Sf, and 6 r ): 

1. scaling of b and /: 

b(r) -> b(r)/5 f , A(r) -> A(r), f(r) -> f(r)/5j, W(r) -> W(r)5 f ; 



18 The "true" horizon position of the physical solution is not r^, but its image under the solution 
generating transformations which we will introduce below. Therefore the value is immaterial 
and can be chosen arbitrarily. 
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2. scaling of b and r: 

b(r) -> 6(r5 6 )/^, A(r) -» A(r5 6 ), /(r) -> /(r<f 6 ), W(r) -»■ W(r5 6 ); 



3. shift in r: 

6(r) -> fo(r - <5 P ), A(r) -> A(r - 5 r ), /(r) -> /(r - 5 r ), W(r) -> W^(r - A 



These transformations leave the system of equations ( |A.1|) invariant, and they 
preserve the regularity condition ( |A.3j ). Thus, they map a regular solution into 
another regular solution. Moreover, they do not act on the function A(r), in particular 
they leave A^ = A(r^) invariant (although they change the other horizon quantities, 
including the horizon position r^). 

The first transformation can be used to set the right normalization of f(r) in the 
UV, f{ruv) — 1- Acting with the second transformation changes the UV behavior of 
the metric and dilaton, or equivalently A (without spoiling the UV value of /); hence 
the strong coupling scale can be fixed by using this symmetry to obtain the same 



UV asymptotic behavior for all solutions, yielding the same value A from eq. ( |4. 11 ) . 



Finally, the shift in r simply serves as a mechanism to fix the same UV boundary 
for all solutions. 

Explicitly, for each A^, given the solution {A(r), b(r), f(r), W(r)}, acting with 
the three transformations above we generate a new solution {A(f), 6(f), f(r), W(r)}: 

~ b (r) = b(r)/(8 b 5 f ) , A(f) = A(r) , (A.4) 
f(r) = f(r)/5j , W(r)=W(r)5 f , 

where f = r5b — S r . For each A^, the parameters 6b(\h), Sf(Xh), S r (Xh) are determined 
such that the UV behavior is the same for all A^. Namely: the UV boundary is at 
f = fuvt the scale factor and dilaton match close to the UV boundary b(fuv + e) = 
buv, X(r uv + e) = \uv (with given constants r uv , buy, Xuv that can be chosen 
arbitrarily 19 and taken to be the same for all solutions). We first set Sf(Xh) such 
that 



then 5 b (X h ) 



and finally 5 r (Xf, 



S f (X h ) = y r f^) , (A.5) 



Sb(X h ) = b{r uv + e)/{b uv 5 f {X h )) (A.6) 



S r (Xh) = r uv 5 b (X h ) - r uv . (A. 7) 



19 Setting buy and \uy determines the energy scale of the theory, namely M, while r uv is 
irrelevant and can be set to zero for convenience. 
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In eg. ( |A.6|) e \ru — ruv\ is chosen such that X{ruv + e ) = Xuv for all Xh, implying 
X{ruv + e) = Xuv with e = e5b <C |fft — fW|. The dependence of the r.h.s of 
eqs.( [A.5| )-( [A.7| ) on Xh is implicit in the dependence of the solutions b(r) = 6(r; X h ), 
X(r) = X(r;X h ),f(r) = f(r;X h ) on X h 

The solutions ( |A.4| ) with the definitions ( A.5| )-( TA~7| ) now represent the physical 
solutions that can be used to compute the thermodynamic properties of the holo- 
graphic plasma, since they all share the same UV behavior: 

r{r uv ) = fuv , f(r uv ) = 1 , (A.8) 
~b(r uv + e) = b uv , X(f uv + e) = X uv , 

for all different values of Xh- 



Next, we describe how to compute the thermodynamic properties of a given 
solution. 

The temperature Inserting into eq.( [2.11| ) the expression for f(rh) = —fhSj 2 
(from now on we drop the ~ to simplify notation), that is obtained after acting with 
the aforementioned transformations on the solutions, we get 



The free energy The free energy is the sum of the two contributions coming from 
Q and C respectively. The entropy contribution C (see eq. ([2.10|) ) is given by the 
area of the black-hole horizon C(Xh) = (b h /5b(Xh)Sf(Xh)) 3 - 

The contribution related to the gluon condensate can be evaluated in two differ- 
ent ways: by directly fitting the metric (or dilaton) fluctuation b a (r) — b(r) close to 
the UV boundary, or by making use of the first law of thermodynamics in the form 
S = —dTj&T. The latter method is more efficient 20 . In particular, summing the 
two contributions in eq.( |2.12| ) yields the following expression for the free energy 

JF{X h )= [ +0a dX h S(X h )^^- . (A.IO) 
h h dA h 



Here S(Xh) is given by (|2.11|) and both the big black-hole and small black-hole 
branches are needed in order the get the full result for the free energy. This is because 
the integral in ([A.IO ) extends to +00, entering deeply in the small black-hole branch. 



In the numeric evaluation the integral is cut off at a large value A^ ^> Xh.min, where 
J 7 (Ah) J^iXh^mm), with Xh,min corresponding to the minimum temperature at 
which the two black-hole branches join. Thus, the value A^ provides a systematic 
way of controlling the numerical error one makes in using the formula ( |A.10| ). 



20 A numerical fit of the fluctuation can be performed in order to compute Q following eq.(pjj). 
It gives larger errors than the integral method described below, though. 
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Thermodynamic functions All thermodynamic functions are given in terms of 
the free energy, entropy and temperature. The energy density e = E/V3 can be read 
from eq. ( |2.8| ), while the pressure p follows from p = —J-/V3. 

The speed of sound is given by (? a = dp/de = s/C v (with C v denoting the specific 
heat at constant volume). 

A.l An alternative technique 

Here we present another numerical method that is used to check the findings ob- 
tained by the previous method above. The main advantage of this method is the 
simplicity. In particular, the idea now is to solve the Einstein's equations using the 
reduced system of variables described in Section 7 of As explained there, the 
information relevant for the physical observables at finite temperature is entirely 
encoded in a coupled system of two first order differential equations for the scalar 
variables X = $'/3A' and Y = f'/AfA'. Note that these functions are invariant un- 
der the diffeomorphisms of the radial coordinate. The idea is to use $ as the radial 
coordinate and solve for X and Y as functions of $. This is enough to construct 
thermodynamical observables. 

Given the potential V^(A) one first integrates the reduced Einstein system, 

from a boundary Ao (chosen to be arbitrarily small) to the horizon A^. In solving 
( |A.ll ) and (|A.12|) one imposes the boundary conditions at the horizon: 



log(A fe /A) 

X^~Y h + 0(\og(\ h /\)), (A.13) 

as A Ah where Yh is given by Yh = 9XhV'(Xh)/32V(Xh). 

Now, the thermodynamic observables can be determined by either of the three 
methods below: 

1. Similar to ( |A.11|) , we solve the same equation that corresponds to zero T that 
is obtained by setting Y = in ( [A.ll ). Let us call the solution of this X . 



One has to specify the boundary condition for numerical integration: Pick an 
arbitrary point A*. Specify the value of X at this point X(A*) = X* such 
that for A ^> A* X asymptotes to the desired value as required from color 
confinement [13]. For the potential fl4.1|) this is X — > — 1/2 for A ^> A*. 
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Having obtained Y, X and Xq one determines the constants of motion Yo and 
Co by fitting the numerical solution to following asymptotic forms: 



Y(A) 
X(X)-X (X) 



Y e ^ (b X) 



-4b 



Yp/2 - Co 

Xa 



+ CqXq 



-16 



(A.14) 
(A.15) 



for A close to the boundary Aq. The constants of motion Co and Y are 
functions of the horizon location A^ (hence temperature) and have the physical 
meaning of energy and enthalpy respectively. To compute them in terms of T 
one uses the relation |T8| , 



T 



Yo 



7T&3(A,. 



-1 4A -> 



(Mc 



-46i 



(A.16) 



Here Aq is the initial value of the scale factor and can be set to an arbitrary 
value as it drops out when we compute thermodynamic observables as functions 
of the ratios T/T c . Similarly £ can be set to 1. 

In summary, one first solves for X, Y and Xq, then determines the constants 
of motion Y and Co from ( |A.14j ) and ( |A.15| ) as functions of A^, next one 
translates them into functions of T using (|A.16| ) and finally obtains the free 
energy density from them using 



A J 



457T 



(6Co(T)-4Y (T)). 



(A.17) 



Here A is the integration constant flj.llQ which should be set by the zero tem- 
perature solution, since we require that the black-hole and the zero T solution 
have the same asymptotics near the boundary. Dependence on A drops out 
when we compute dimensionless quantities such as s/T 3 . T c is obtained by the 
value of T that makes ( |A.17 ) zero. Given ( |A.17| ) all the rest of the thermody- 



namics follow from the formulae presented in section 7 of |T8|. 



There is a second method that is simpler. Here, one does not need to determine 
Xq. One solves for X and Y as above, than uses the first law of thermodynamics 
to construct T directly from the entropy S, using ( |A.10| ). Here S in ( |A.10| ) 
is given by the Bekenstein-Hawking formula ( |2.11| ) where the scale factor is 
constructed from X by (T8| : 



A{X) = A + 



Ao 



1 dX 



(A.18) 



whereas the T in ( |A.10|) is given by (|A.16|) above. T c is determined by the 
value that makes ( |A.10| ) zero. 
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3. We end this subsection by presenting an even simpler formula, that side-steps 
most of the complications above and determines the thermodynamic functions 
semi- directly, in terms of the potential. Given V, one solves for X as above. 
Then the entropy density s = S/V 3 N^ is given by 



s 

j>3 



768tt 2 



e -4/^fX(A) 



5£ 6 



V(X, 



16384tt 2 f W{\ h ) 
1215£3 \V{\) 



(A.19) 



In the second line we used the definition of the thermal superpotential in terms 
of X, \TB\. We refer to |TH| for a derivation of ( |A.19 ). Using this and ( A. 16 ) one 
easily obtains the thermodynamic functions. We note that the dimensionless 
quantity s/T 3 is one of the main functions that is compared with the lattice 
data, hence this formula provides a semi-direct method of relating the data to 
be compared with lattice and the dilaton potential. The reason for its being 
semi-direct rather than direct is because one still has to determine T c in order 
to set the scale which requires use of (|A.19|) in (|A.10|). 
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